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Abstract 

The ultra-high energy cosmic rays recently detected by several air shower experiments 
could have an extragalactic origin. In this case, the nearest active galaxy Centaurus A 
might be the source of the most energetic particles ever detected on Earth. We have used 
recent radio observations in order to estimate the arrival energy of the protons accelerated 
by strong shock fronts in the outer parts of this southern radio source. We expect detections 
corresponding to particles with energies up to ^ 2.2 x 10^^ eV and an arrival direction 
of (/ « 310°, b « 20°) in galactic coordinates. The future Southern Hemisphere Pierre 
Auger Observatory might provide a decisive test for extragalactic models of the origin of the 
ultra-high energy cosmic rays. 
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1 Introduction 



The origin of cosmic rays (CR) detected in the neighborhood of Earth is still an unsettled issue. 
Most of the low energy particles are believed to be produced by supernova explosions in the 
galactic interstellar medium Higher energy CR could be accelerated by stellar explosions in 
strong stellar winds H]. The basic acceleration mechanism in all these models is diffuse shock 
acceleration Q. Different observational facts, such as a break in the CR spectrum at ~ 5 x 10^^ 
eV, the change in the composition of the rays from heavy nuclei below this break to light nuclei 
and protons at higher energies, and the arrival directions of the highest energetic particles seem 
to support the idea that the component with energies below the break is predominantly galactic 
whilst the particles with energies above 10^^ eV have an extragalactic origin. The hot spots 
of extended radio galaxies (i.e. the end-points of powerful jets ejected by the central active 
galactic nucleus) have been suggested as the dominant source of ultra-high energy CR Q. In 
such a case, particles are assumed to be accelerated by first order Fermi mechanism through the 
strong, nonrelativistic shock waves generated by huge plasma collisions. 

When the ultra-high energy protons are injected in the intergalactic medium, traveling losses 
originated by interactions with the 2.73 K cosmic background radiation are expected to produce 
a high energy cutoff in the particle spectrum at energies lO'^*' eV, according to the predictions 
of Greisen and Zatsepin and Kuzmin Q (the so-called "GZK cutoff"). Nevertheless, recent 
detections with the Akeno Giant Air Shower Array and the Utah Fly's Eye detector |^, ^ 
have revealed the existence of CR with energies above the mentioned cutoff. A detailed analysis 
of the arrival directions of these and other energetic events has shown a correlation with the 
general direction of the supergalactic plane, where many potential extragalactic sources of CR 



are located [10|. Despite these correlations, there seems to be no conclusive evidence for an 

extragalactic origin of the most energetic CR. 

In this paper we present a concrete candidate for testing the acceleration models that involve 

hot spots and strong shock waves as the source of the ultra-high energy component of the 

CR. Using polarization measurements and synchrotron emission as a trace, we can make an 

estimate of the maximum energy of the protons injected in the extragalactic medium and the 

corresponding arrival energy for the case of the nearest active galaxy: the southern radio source 

Centaurus A (Gen A)Q. If the basic model is correct, a southern hemisphere detector should 
^Cen A was suggested as a possible source of ultra-high energy cosmic rays by G. Cavallo from quite 
general energetic arguments. 
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register events produced by particles with energies beyond lO^i eV as we shall see. But first, let 
us briefly discuss some physical features of Cen A. 

2 Cen A and CR acceleration 

Cen A is a complex and extremely powerful radio source identified at optical frequencies with 
the galaxy NGC 5128. It is the nearest radio galaxy, at a distance of ^ 3.5 Mpc according to 



recent works on globular clusters |T3] and planetary nebulae |13|. Radio observations at different 



wavelengths [14, 15 1 show a structure composed by a compact core, a one-sided jet. Double Inner 
Lobes, a Northern Middle Lobe, and two Giant Outer Lobes. This morphology, together with 
the polarization data obtained by Junkes et al. [^] and the large-scale radio spectral index 
distribution computed by two of us |jl^, strongly support the picture of an active radio galaxy 
with a jet forming a relatively small angle with the line of sight. The jet would be responsible 
for the formation of the Northern Inner and Middle Lobes when interacting with the interstellar 
and intergalactic medium, respectively. The Northern Middle Lobe can be interpreted as a 



"working surface" |17] at the end of the jet, a place where strong shocks are produced by plasma 
collisions, i.e. it can be considered as the hot spot of a galaxy with a peculiar orientation. 

The acceleration of particles in the hot spot is the result of repeated scattering back and 
forth across a strong shock front in a partially turbulent magnetic field. This process has been 



studied in detail by Biermann and Strittmatter [18|. Assuming that the energy density per unit 



of wave number of MHD turbulence is of Kolmogorov type, i.e. I{k) a k ^ with s = 5/3, we 
have an acceleration time scale for protons given by: 

where is the jet velocity in units of c, u is the ratio of turbulent to ambient magnetic 
energy density in the hot spot (of radius R^^), and B is the total magnetic field strength. The 
acceleration process will be efficient as long as the energy losses by synchrotron radiation and 
photon-proton interactions do not become dominant. Considering an average cross section a-yp 
for the three dominant pion-producing interactions |19]: 

7 + p ^ p + 7r° 
7 + p — > n + vr"'" 
7 + p^p + 7r++7r~, 



3 



we get the time scale of the energy losses within a certainty of 80%: 

dt J^^^^ ^ a^mjB^l + aA) P ^ ' 

where a stands for the ratio of photon to magnetic energy densities, cij, is the classical Thom- 
son cross section, and A gives a measure of the relative strength of 'jp interactions against the 
synchrotron emission. Biermann and Strittmatter ||l^ have estimated A ~ 200, almost inde- 
pendently of the source parameters. The most energetic protons injected in the intergalactic 
medium will have an energy that can be obtained by balancing the energy gains and losses: 

^p,max = 1259.3 C f^Y' ^ R^-l/^ (1 + ^«)-3/4 (3^ 

V / me 

The acceleration model also predicts a power-law spectrum for the number of particles per unit 
of energy: N{E) oc E'P, {E < Ep,^^,). 

In the case of Cen A, we can estimate u from the radio spectral index of the synchrotron 
emission in the Northern Middle Lobe and the observed degree of linear polarization in the same 
region. We get u ~ 0.4. The size of the hot spot can be directly measured from the large-scale 
map obtained by Junkes et al. with the assumed distance of 3.5 Mpc, giving as a result 
i?jjg f« 1.75 kpc. 

The equipartition magnetic field can be obtained for the pre-shock region from the detailed 
radio observations by Burns et al [^]. The field component parallel to the shock will be amplified 
in the post-shock region by a compression factor ^. In the case of a strong, nonrelativistic shock 
front, ^ ^ 4 [^], and then, if ~ we have: 

B ^ {e'^ + ifi^B^ 

« 5 X 10~5 G 

An enhancement of the B\\ component in the Northern Middle Lobe can be clearly seen in 
the polarization maps displayed in ref. [15|. The value of /^j^x is uncertain. We shall assume 



s::: 0.3, a typical value for similar sources |21|. 



With the above mentioned values for the input parameters in Eq. (|^, the maximum energy 
of the protons injected in the intergalactic space results 

£^p,max~2.7x lO^i eV (4) 

We can infer the index in the power-law spectrum from multifrequency observations of the 
synchroton radiation produced by the leptonic component of the particles accelerated in the 



source (see the standard formulae, for instance, in the book by Pacholczyk [^]). Using the 
radio spectral index obtained by Combi and Romero for the hot spot region [^], we get p = 2.2. 

3 Propagation effects 

In this section we shall briefly discuss the proton energy losses due to collisions with microwave 
photons during their travel to Earth, and the modifications in the energy spectrum arised in 
these interactions. The losses can be estimated as: 



1 — exp 



2-fpkT 



de (5) 



where T=2.73 K, r is the traveling time, fjie) is the cross section weighted by inelasticity of the 
jth reaction channel, and the sum is carried out over all channels. The integration is over the 
photon energy in the proton rest frame, from the threshold. For energies higher than 3 x 10^'^ eV, 
and taking into account just the dominant pion-producing interactions, this equation becomes 

m-- 

-^^^1.8x10-^ yr-i (6) 
Ep dr 

According to this, the losses are almost neglectable due to the proximity of the source, yielding 
a maximum arrival energy 

^p'max ~ 2.2 X 10^1 eV . The arrival direction of these CR will be 
approximately (/ ~ 310°, b ~ 20°), in galactic coordinates. Intergalactic and galactic magnetic 
fields with strength of the order of nG and fxG, respectively, cause no relevant deviations from 
the line of sight, due to the distances and energies involved in the case under consideration. 
However, a spreading in the arrival direction pattern should be observed for a relatively large 
number of low energy events coming from Cen A. 

Both the JTT energy losses and the conservation of the total number of protons in the spec- 
trum are expected to produce only a slight modification of the original energy spectrum. This 
modification can be quantitatively described in terms of a factor given by 

In the case of a nearby source, an analitical expression of i]{Ep, r) has been found by Berezinsky 
and Grigor'eva |23|: 
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At very high energies {Ep > 5x10^*^), d\nP{Ep)/d\nEp 0, and then r]{Ep,T) 1 — {p—1)tujq, 
where u>o ~ 1.8 x 10^^ yr^^. Since rj < 1, there wih be a shght fall in the spectrum, as we already 
mentioned in the case of the highest energy protons. At lower energies (Ep^Sx 1020 eV) the 
modification factor is 

where ejr ~ 3 x 10^^ gy^ and t-j^ ~ 1.7r. This will produce a small bump in the spectrum. In 
brief, the overall spectrum will be scarcely modified. 

An additional interesting consequence of the model is that the 7 + p — > p + vr^ interactions 
would be a source of 7-rays through decay (tt^ —5-27). Considering a power-law dependence 
for the isotropic and homogeneous photon distribution in the source [|l^]: 

N^{e) = no (j-^ eo<e< e^^ax (10) 



we can compute the energy range of the total pion spectrum |24]. Using 7p,max ~ 10^^ and 
Cmax ~ 1 MeV, we obtain (in the observer's frame): 

E^e[4x 10^°eV, 5 x lO^^eV] (11) 

If the proton distribution N(Ep) has the same index as the electron distribution N(Ee), as we 
have assumed, the pion distribution function will be 

N{E^) oc E'"^ (12) 

In the case of Cen A, we have N^E^^) oc {'^^^ m-,^ c^)"-"^-^. 

Emission at low 7-ray frequencies from Cen A has been detected by the Compton Gamma 
Ray Observatory ||2^, and by experiments with balloons |^^. However, these 7-rays could not 
be produced by vr" decays in the acceleration region because, according to (11), the relevant 
energy ranges from ^ 20 GeV up to a value comparable to those of ultra-high energy CR. It 
might be worthwhile then to perform observations of Cen A with instruments sensitive to TeV 
energies, for instance with the new Imaging Atmospheric Cerenkov Telescopes |27]. 



4 Final comment 

We expect the detection of an extensive air shower produced by the collision of protons with 
energies Ep > 100 EeV with the Earth's atmosphere. Since the arrival direction of these ex- 
tremely energetic particles is well known, the future Hybrid Pierre Auger Observatory ||2^, ^ 
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will provide a decisive test for the models that invoke an extragalactic origin for the ultra-high 
energy component of the cosmic rays. 
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